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On approaching the glass transition, the microscopic kinetic unit spends increasing time rattling in the cage 
of the first neighbours whereas its average escape time, the structural relaxation time Tq, increases from a 
few picoseconds up to thousands of seconds. A thorough study of the correlation between Tc and the rattling 
amplitude, expressed by the Debye-Waller factor (DW), was carried out. Molecular-dynamics (MD) simulations 
of both a model polymer system and a binary mixture were performed by varying the temperature, the density 
p, the potential and the polymer length to consider the structural relaxation as well as both the rotational and 
the translation diffusion. The present simulations, together with MD studies on other glassformers, evidence 
the scaling between the structural relaxation and the caged dynamics. An analytic model of the master curve is 

developed in terms of two characteristic length scales ^ and a^2 pertaining to the distance to be covered 
by the kinetic unit to reach a transition state. The model does not imply Ta divergences. The comparison with 
the experiments supports the numerical evidence over a range of relaxation times as wide as about eighteen 
orders of magnitude. A comparison with other scaling and correlation procedures is presented. In particular, 

the density scaling of the length scales a^^^^, oc p~^^^ is shown to be not supported by the present 

simulations. The study suggests that the equilibrium and the moderately supercooled states of the glassformers 
possess key information on the huge slowing-down of their relaxation close to the glass transition. The latter, 
according to the present simulations, exhibits features consistent with the Lindemann melting criterion and the 
free-volume model. 

PACS numbers: 64.70.Q-,02.70.Ns 
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I. INTRODUCTION 

When they are cooled or compressed, several systems like 
liquids, mixtures, polymers, bio-materials, metals and molten 
salts may avoid the crystallization and, following a huge in- 
crease of the viscosity, finally freeze into a glass, a microscop- 
ically disordered solid-like state. Understanding the extraor- 
dinary viscous slow-down that accompanies glass formation 
is a major scientific challenge OlUjllt]. 

On approaching the glass transition (GT), trapping effects 
are more and more prominent. The average escape time from 
the cage of the first neighbors, i.e. the structural relaxation 
time Tq, increases from a few picoseconds up to thousands 
of seconds. The rattling motion inside the cage occurs on pi- 
cosecond time scales with amplitude (u^)^/^, the so called 
Debye-Waller factor (DW). The DW factor is clearly related 
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to the short-time elastic properties of the systems At 
first sight, due to the extreme time-scale separation between 
the rattling motion (~ 10"^^s) and the relaxation (~ VS^s 
at GT), one expects the complete independence of the two 
motions. However, already in 1943 Tobolsky, Powell, and 
Eyring pointed out that there could be a relation between the 
curvature of the potential well near the minimum (controlling 
the DW factor) and the height of the energy barrier (limiting 
the flow process), thus establishing a relation between the in- 
stantaneous shear modulus Goo and the shear viscosity 77 isj] . 
Later, the diffusive motion was described as natural conse- 
quence of the dynamic equilibrium between vibrational and 
configurational quantum states |6] and the free-energy barrier 
for viscous flow was found as being proportional to Goo (T) 

A firmer basis to connect fast and slow degrees of free- 
dom was developed by Hall and Wolynes who, assuming that 
atomic motion is restricted to cells, pictured the glass transi- 
tion as a freezing in an aperiodic crystal structure (ACS) mod- 
eled by the density functional theory |8]. As a result, the vis- 
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cous flow is described in terms of activated jumps over energy 
barriers AE oc kgTa'^/ {v?) where a is the displacement to 
reach the transition state and ks the Bohzmann constant. The 
usual rate theory leads to the Hall-Wolynes equation (HW): 



(1) 



The ACS model is expected to fail when Tq. becomes com- 
parable to the typical rattling times of each atom in the 
cage of the surrounding atoms, corresponding to picosecond 
timescales. That condition is quite mild, e.g. in Selenium it 
occurs at + IQAK (Tm is the melting temperature) J^. 
Buchenau and Zorn derived a relation very similar to EqlUin 
terms of soft vibrational modes |9]: 



Tq, , oc exp 
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where uq is a critical displacement to allow for the elemen- 
tary flow or a-relaxation process and {v?)ioc is the difference 
between the DW factor in the liquid phase (u^) and its ex- 
trapolation from the low-temperature values. The definition 
of (u^)/oc affects the plot log 77 vs. 1/ {u^)ioc- If the extrap- 
olation of either the glass or the crystal contribution is sub- 
tracted from the DW factor of selenium, a convex curve or 
a straight line are seen, respectively [9]. The fact that many 
glass-formers have no underlying crystalline phases, as well 
as the fact that in other studies removing the glass contribu- 
tion, differently from selenium, the plot log rj vs. l/{u'^)ioc 
is a straight Une ifTol [Till , raises some ambiguities about the 
above subtractions. Buchenau and Zorn also noted that, if no 
subtraction is made, the curve log 77 vs. l/{u^) for selenium 
is concave, namely the HW equation, EqlT] is not obeyed. 
The HW equation has been derived in the framework of the so 
called elastic models (for a review see ref.il), like the shov- 
ing model Ijj . 

The HW equation states that the glass softens when the DW 
factor exceeds a critical value, which is reminiscent of the 
Lindemann melting criterion for crystalline solids |14], The 
empirical law Tg ~ 2/3r„j (Tg is the glass-transition temper- 
ature) id i [H also suggests that the melting and the glass 
transition have a common basis. This viewpoint led to an al- 
ternative derivation of EqlT] 1 16] and motivated extensions of 
the Lindemann criterion to glasses IITtIi . The closeness of the 
HW equation with free-volume concepts ifTsIl was noted 
and investigated numerically [ 19]. 

Other studies noted a relation between the fast vibrational 
dynamics and the long-time relaxation both far ll20l Ull] and 
close to the glass transition ||3 [isl 113, HH IH. A numeri- 
cal investigations pointed out that the short-time DW hetero- 
geneities predict the spatial distribution of the long-time dy- 
namic propensities [20]. The fragility tti, a steepness index 
of how fast the viscosity rj or Tq, increase close to Tg [1], has 
been also considered. It has been proposed that variations in 
the fragility originate in differences between their vibrational 
heat capacities, harmonic and anharmonic [22] and depend 
on changes in the vibrational properties of individual energy 
minima of the energy landscape in addition to their total num- 
ber and spread in energy [21]. The temperature dependence 



of the DW factor around the glass transition has been also 
studied ^ [Tsl [H. It was seen that for strong glass- 
formers (small fragility) DW is almost linear with tempera- 
ture, whereas a stronger than linear dependence takes place 
for fragile systems pointing to increasing anharmonicity of the 
short-time dynamics. With a distinct approach further studies 
established correlations between the vibrational dynamics and 
the relaxation close to the glass transition, as quantified by the 
fragility iH [M IH S with controversies jH. Finally, 
as further examples of studies comparing the fast and the slow 
dynamics, we point out the correlations between the structural 
and the secondary relaxations in a supercooled liquid [31], as 
well as between the apparent activation energy above the glass 
transition and the fragility ll28ll29ll . 

In a recent paper we reported the universal dependence be- 
tween the structural relaxation time and the DW factor for 
a model polymer ll32ll . The universal scaling curve, which 
is described by a simple generalization of the HW equation 
(EqUll, fits with the existing experimental data from super- 
cooled liquids, polymers and metallic glasses over about eigh- 
teen decades of relaxation times and a very wide range of 
fragilities. Here we show by novel numerical simulations that 
the scaling holds for binary mixtures with different interacting 
potentials, density and temperatures, i.e. for an atomic, het- 
erogeneous system different by the molecular, homogeneous 
one considered in ref.[32[. Moreover, we prove that it holds 
not only for the translational degrees of freedom but for the 
rotational ones of the model polymer as well. Comparisons 
with other numerical studies if ssl [34l Issi [36l IjtIi . novel ex- 
perimental data [38] as well as other scaling and correlation 
procedures is also presented. 

The paper is organized as follows: in Sec|ll]the HW equa- 
tion is suitably generalized. In Sec|III]the numerical meth- 
ods are described. The results are presented and discussed in 
SecjIVI The conclusions are given in SecfVl 



II. GENERALIZED HALL-WOLYNES EQUATION 

One basic assumption of the original HW equation, EqlT] is 
that the distance to reach the transition state has a characteris- 
tic value a. Actually, this length scale is dispersed. To model 
the related distribution, it is assumed that the latter does not 
depend on the state parameters such as the temperature, the 
density or the interacting potential. This complies with the 
spirit of ref.[!8|| where the a distance is said to be mostly con- 
trolled by the geometrical packings. It is also known that, irre- 
spective of the relaxation time Tq, the average distance moved 
by the relaxing unit within is about the same, i.e. a fraction 
of the molecular diameter [j]]|. As a suitable choice, the distri- 
bution of the squared distances p{o?) is taken as a truncated 
gaussian form 




(a"-a^)- 
2a\ 



if ^ CLfnin 

Otherwise 
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where A is the normalization and a^^^^ is the minimum dis- 
placement to reach the transition state. Averaging the HW 
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EqlUover the distribution given by Eq|3] yields the following 
generalized HW equation (GHW): 



= B7V[(it^)]exp 



2(u2) 8(u2)2 



(4) 

(5) 



where i? is a constant and the normalization factor M[{u'^)] 
reads: 
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If a2 > al^i^, 1 < N[{u^)] < 2, namely 7V[(u2)] depends 
very weakly on the DW factor (u^) and the influence of the 
truncation is negligible. Then, tq = BAf[{u^)] ~ const and 
Eq|5]reduces to: 



Ta = To exp 



(7) 



An analogous law holds for the viscosity rj. Owing to the finite 
value of the DW factor, Eq|7]does not imply the divergence of 

Ta MM- 

The motivations behind the gaussian form of p{a'^) mainly 
rely on the Central Limit Theorem. In fact, (r§ in the no- 
tation of ref.|8]) is the cumulative displacement of the Nm 
particle that move |8]. Other supporting facts for the gaus- 
sian form of p(a2) are the following. If the kinetic unit per- 
forms harmonic oscillations around the equilibrium position 
with an effective spring constant k, the DW factor becomes 



ksT/k and Eq. dTji reduces to: 
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The above expression was reported for both supercooled liq- 
uids ll4lll and polymers |42]. Along a similar line of reason- 
ing, assuming harmonic oscillations leads to the following ex- 
pression for the energy barrier height AE 1^: 



AE = -ko? 



(9) 



Eq|9]allows one to reinterpret the gaussian form of p{a?) as a 
gaussian distribution of energy barriers [43]. Substituting Eq. 
(|9]l into Eq. ([8]), one recovers a key result of the facilitated 
model of glass-formers developed by Garrahan and Chandler 
iH: 



To exp 
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m. METHODS 
A. Models 

1. Polymer melt 

A coarse-grained model of a linear polymer chain is used. 
Torsional potentials are neglected. We considered a system 
of N,n — 2000 monomers in all cases but AI — 3 where 
Njn — 2001. Non-bonded monomers at a distance r interact 
via the truncated parametric potential: 



p- q 



P 



q 



Uc 
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where and the value of the constant Ucut is cho- 

sen to ensure Up^r) = at r > Tc = 2.5(T. The minimum 
of the potential Up^r) is at r ~ a* , with a constant depth 



e. Note that Uq^r) = Up^r). Bonded 



monomers interact with a potential which is the sum of the 
FENE (Finitely Extendible Nonlinear Elastic) potential and 
the Lennard-Jones (LJ) potential |45]. The resulting bond 
length is 6 = 0.97(t within few percent. We set cr = 1, e = 1. 
The time unit is tmd = (jna^ j e)^!'^ , with m being the mass 
of the monomer Temperature is in units of e/fc^, where ks 
is the Boltzmann constant. We set m = fc^ = 1. NPT and 
NTV ensembles have been used for equilibration runs while 
NVE ensemble has been used for production runs for a given 
state point (labelled by the multiplets {T, p, Af,p, q\). NPT 
and NTV ensembles have been simulated with the extended 
system method introduced by Andersen |46] and Nose |47l]. 
The numerical integration of the augmented Hamiltonian has 
been performed through the reversible multiple time steps al- 
gorithm (r-RESPA algorithm) et a/. |48]. In particular, the 
NPT and NTV Liouville operators have been factorized us- 
ing the Trotter theorem [49.1 separating the short range and 
long range contributions of the potential Up^r) (see Eq. fTTT l. 
according to the WCA decomposition llsoll . 



2. Binary mixtures 

An 80:20 binary mixture (BM) of Nbm = 1000 particles is 
considered. The two species are labelled A, B and particles 
interact via the potential: 



(r) 



p-q 



P 



'a,p 



(12) 

that is similar to Eq. ( fTTT l. except that the well height and 
the minimum of the potential now depend on the interacting 
species, being a, (3 £ A, B with aAA = 1-0, aAB = 0.8, 
(^BB = 0.88, eAA = 1-0, eAB = 1-5, ess = 0.5. Note that 
setting q = 12, p = 6 in Eq. ( fT2] i and the above choices for 
aa.f3 and ea.fs, reduce the model to the well-known LJ Kob- 
Andersen model (BMLJ)|51, 52, 53]. The system was equi- 
librated in the NTV ensemble and the production runs were 
carried out in the NVE ensemble. NTV runs used a standard 
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Nose method (4^. The "velocity verlet" integration algorithm 
was used both in the NVE and NVT ensembles I.54i1 . 



IV. RESULTS AND DISCUSSION 



A. Relaxation and transport properties 



First, the monomer dynamics has been studied. To this aim, 
one defines the mean square displacement (MSD) (r^(t)) as: 



l^(||x.(t)-x.(0)|p) 



(13) 



In addition to MSD the self part of the intermediate scattering 
function (ISF) is also considered: 



N 



1 



(14) 



ISF was evaluated at q = Qmax, the maximum of the static 
structure factor. N — and N = TV^m for polymers and 
binary mixtures, respectively, is the position of the i-th 
monomer (polymers) or particle (BM). Note that for BM MSD 
is averaged over both species A and B. 

FiglU shows typical MSD and ISF curves of the polymeric 
monomers. At very short times (ballistic regime) MSD in- 
creases according to = (iksT /m)t^ and ISF starts 
to decay. The repeated collisions with the other monomers 
slow the displacement of the tagged one, as evinced by the 
knee of MSD at < \/T2/r2o ^ 0.17, where fio is an effec- 
tive collision frequency, i.e. it is the mean small-oscillation 
frequency of the monomer in the potential well produced by 
the surrounding ones kept at their equilibrium positions ifssll . 
At later times a quasi-plateau region, also found in ISF, oc- 
curs when the temperature is lowered and/or the density in- 
creased. This signals the increased caging of the particle. The 
latter is released after an average time t^, defined by the re- 
lation Fs{qmax, Ta) = 6^^. For t > Ta MSD increases more 
steeply. The monomers of short chains (M < 3) undergo dif- 
fusive motion {r^{t)) cx with 5=1. For longer chains, 
owing to the increased connectivity, the onset of the diffusion 
is pr eceded by a subdiffusive region (S < 1, Rouse regime) 

The monomer dynamics depends in a complex way on the 
state parameters. Nonetheless, if two states (labelled by mul- 
tiplets {T, p, M, p, q}) have equal relaxation time r^, the cor- 
responding MSD and ISF curves coincide from times fairly 
longer than down to the crossover to the ballistic regime 
and even at shorter times if the states have equal temperatures. 
Examples are shown in FigH] See EPAPS supplementary ma- 
terial at [URL will be inserted by AIP] for the details on all 
the investigated states. Notice that the coincidence of MSD 
and ISF curves of states with equal at intermediate times 
(t ^ Ta) must not be confused with the customary superpo- 
sition of ISF curves at long times (t > r^) following a suit- 
able logaritmic time shift (see the lower-panel inset of Fig[T]). 
States with coinciding MSD and ISF have close non-gaussian 
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FIG. 1: Monomer dynamics in the polymer melt. Top: MSD time- 
dependence for polymers in selected cases. MSDs are multiplied by 
indicated factors. Inset: corresponding MSD slope A{t); the uncer- 
tainty range on the position of the minimum at t* = 1.0(4) (black 
line) is bounded by the vertical colored lines. Bottom: correspond- 
ing ISF curves for polymers. Inset: superposition of the ISF curves. 
Four sets of clustered curves (A through D) show that, if states have 
equal Tq (marked with dots on each curve), the MSD and ISF curves 
coincide from times fairly longer than Ta down to the crossover to 
the ballistic regime at least. 



. This is shown by the non-gaussian parameter 



properties (S; 
(NGP): 



3 jrHt)) 



(15) 



where (r'*(t)) is defined analogously to MSD, Eq[T3] Plots of 
a2{t) for the same states of Fig[T]are shown in Fig|2l One no- 
tices that states with coinciding MSD and ISF have coinciding 
NGP as well. 

Owing to the fully-flexible character of the chain, the struc- 
tural relaxation time little depends on the chain length M iIstIi . 
Much stronger dependence is expected for the diffusion coef- 
ficient and the reorientation time of the whole chain which for 
unentangled chains (A/ < 32 [45]) scale as D^^ cx AI and 
Tee cx; A/^, respectively Isstl . These processes set the long- 
time dynamics of the chain and it is interesting to see if states 
with coinciding MSD, ISF and NGP, involving short and in- 
termediate time scales, also exhibit coinciding translational 
and rotational diffusion. To this aim, the global rotational dy- 
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FIG. 2: The non-gaussian parameter of the polymer states plotted in 
FigE 

namics of the chain has been investigated by the correlation 
function of the end-to-end vector: 

Ceeit) = -T— ^ W • I^^W) (16) 

where Ri(i) is the vector joining the first and the last 
monomer (i.e. the end monomers) of the i-th polymer in the 
melt and 

2—1 

Cee{t) monitors the collective relaxation, whereas ISF the 
single-particle one. One defines T^e via the equation 
Cee(Tee) = 1. For uncntanglcd polymers T^e ^ AAPt^ 
llssll . FiglDplots the correlation function Cee in dependence 
of the reduced time At/A-'P, i.e. the curves are scaled onto 
the one of the dimer (M = 2) whose end-to-end vector is 
the bond itself. Having removed the chain-length dependence 
by proper rescaling, the states with coinciding MSD, ISF and 
NGP exhibit coinciding end-to-end correlation loss Cee{t) 
too. Note that the polymer states contributing to one clus- 
ter of scaled curves have not necessarily equal chain length 
(see EPAPS supplementary material at [URL will be inserted 
by AIP] for the details on all the investigated states). This 
is also evidenced by inspecting FiglTftop). In fact, up to 
t ^ Ta the connectivity effects are negligible and, irrespec- 
tive of the M value, MSD curves coincide, whereas at longer 
times monomer bonding comes into play and the curves start 
to differ from each other due to the different chain lengths 
ifssll . It must be noted that, since D oc R^e/Tee Ull, the col- 
lapse of the correlation function Cee ensures that the quantity 
D ■ M is identical for states with coinciding MSD, ISF and 
NGP 

We now consider the BM system. Fig|4]shows typical MSD 
and ISF curves. Note that for the BM system these quan- 
tities are averaged over both A and B species. One sees 
that, if two states (labelled by multiplets {T, p,p,q}) have 




FIG. 3: Correlation functions of the end-to-end vector of the states of 
FigH The scaled time removes the chain length dependence. Poly- 
mer states contributing to one cluster of scaled curves have not nec- 
essarily equal chain length. Dots mark the time iree/M'^. 



equal relaxation time Tq, the corresponding MSD and ISF 
curves coincide at least from the end of the ballistic regime 
onwards, i.e. the states have equal diffusion coefficients 
D ~ limt^oo{r'^{t))/6- This shows that, due to the missing 
connectivity of BM, the MSD coincidence is not interrupted 
at < > Tq as it happens in polymers (see FiglT]l. Remarkably, 
if the temperatures of the BM states are the same the MSD 
and ISF coincidence include the ballistic regime too. In full 
analogy with the polymer case Fig|5] shows that states with 
coinciding MSD and ISF have coinciding NGP too. 



B. Scaling between relaxation and caged dynamics 

1. Polymers and Binary Mixtures 

The results of Sec lIV Al stronglv support the conclusion of 
a close correlation between the caged dynamics at short times 
and the long-time dynamics, including both the structural re- 
laxation, the chain reorientation and the diffusivity. In order 
to better evidence such a correlation a suitable metrics of the 
caged dynamics is needed. This is achieved by considering the 
Debye-Waller (DW) factor (u^), a characteristic length scale 
of the particle temporarily trapped into the cage. 

Preliminarily, one has to clarify if the cage exists. From this 
respect, it must be pointed out that the product JIoTq is ^ 20 
for states with the fastest structural relaxation, meaning that 
the structural relaxation time is at least one order of magnitude 
longer than the collision time. Furthermore, in the present 
study the time velocity correlation function (VCF), after a first 
large drop due to pair collisions, reverses the sign since the 
monomer rebounds from the cage wall (data not shown). 

The DW factor is a characteristic length scale of the rattling 
motion into the cage. The measure of the DW factor must 
take place in a time window where both the inertial and the 
relaxation effects are not present. To clearly identify that time 
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FIG. 4: Average particle dynamics in the BM system. Top: MSD 
of selected cases. MSDs are multiplied by the indicated factors. The 
dashed line shows the position of the minimum of A{t), which is 
plotted in the inset. Note that it shifts at longer times for states 
with slower relaxation. Bottom: corresponding ISF curves. Inset: 
superposition of the ISF curves. Four sets of clustered curves (A 
through D) show that, if states have equal Ta (marked with dots on 
each curve), MSD and ISF curves coincide at least from the end of 
the ballistic regime onwards. If the temperatures are the same, the 
coincidence include the ballistic regime too. 



window we consider the slope of MSD in the log-log plot 

d\og{r^{t)) 



A{t) 



dlogt 



(18) 



Representative plots of A(<) for the polymer system are given 
in the top inset of FiglUand Fig|4] A{t) exhibits a clear min- 
imum at <* = 1.0(4) (corresponding to an inflection point in 
the log-log plot of (r^(i))) that separates two regimes. The 
short- and the long-time limits of A(t) correspond to the bal- 
listic (A(0) — 2) and the diffusive regimes (A(oo) — 1), 
respectively. At short times, t < 0.7 < t* the inertial effects 
become apparent. At long times (i > > t*) relaxation 
sets in. It may be shown that a minimum of A(t) implies that 
VCF exhibits a negative tail at long times. A monotonically 
decreasing VCF, i.e. with no cage effect, leads to a mono- 
tonically decreasing A{t). Therefore, MSD at t* is a mean 
localization length and the DW factor is defined as 



FIG. 5: The non-gaussian parameter of the BM states of Fig|4l 



Notice that t*, corresponding to about 1 — lOps 116011 . is con- 
sistent with the time scales of the experimental measurement 
of the DW factor, e.g. see l'9'l. As far as BMs are concerned, 
the top inset of Fig|4]shows that the time dependence of A{t) 
is rather similar to the polymer case. Fig|6](top) plots the po- 
sition of the minimum of A{t), t*, for the polymer and BM 
systems. The plot shows that it is virtually constant in poly- 
mers, whereas it increases with the structural relaxation time 
in BM (see also Fig|4]top panel). 

The DW factor is usually experimentally measured by using 
ISF and considering the height h of the plateau signalling the 
cage effects (see FiglUand Fig|2]l via the relation: 



(20) 



(19) 



where h is the ISF height at the inflection point of the plateau. 
FiglSbottom) shows that (u^)/sf and {v?) msd = {u^) are 
quite close to each other with constant ratio within our ac- 
curacy. This will have important consequences when com- 
paring the MD simulations with the experimental results in 
Sec ]IVB3l 

To make it explicit the correlation between the relaxation 
and the caged dynamics, Fig|2]shows the dependence of both 
Tq and the scaled average chain reorientation time T^e on the 
DW factor The data collapse on two well-defined master 
curves. The one concerning is well fitted by Eq|2] The 
master curve of the scaled r^e is different. At large DW factor, 
i.e. fast relaxation, the ratio AT^e/ (ApTa) is roughly constant 
and decreases when the relaxation slow down, as previously 
reported ll6Tll . The behaviour at large DW factor is consistent 
with the Rouse theory concerning the dynamics of unentan- 
gled polymers stating that the different relaxation time scales 
are proportional to each other fssll . The Rouse theory also pre- 
dicts the scaling r^e oc AP which is indeed observed even for 
states with very slowed-down dynamics, see Figl3] However, 
the differences between the master curves for the chain re- 
orientation and the structural relaxation, which become more 
apparent for states with sluggish relaxation, evidence one ba- 
sic limit of that theory, i.e. the assumed gaussian and homo- 
geneous character of the monomer displacements. One an- 
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FIG. 7: The structural relaxation time Tq of polymers and BM and 
the scaled average polymer reorientation time Tee vs the DW factor 
{u^). Empty circles highlight the cases plotted in Fig[T] For clarity 
sake the BM cases plotted in Fig[2]are not pointed up. The dashed 
line across the Tq curve is Eq. |2]log = a + j3{u^)~^ + 'y(u^)~^ 



with a 



-0.424(1),/? = aV(21nlO) = 2.7(1) • 10" 



,7 = 



CTf2/(81nlO) = 3.41(3) ■ lO'-^*. Additional data on the collective 
relaxation time r are also plotted (V) 1 19]. The dotted curve is ob- 
tained by vertically shifting the dashed curve (a' = a + 0.205(5)). 
The dashed curve across the scaled average chain reorientation time 
curve is a guide for the eyes. Inset: the maximum of the non-gaussian 
parameter a2 max of the A-F clusters of polymer states vs. the ratio 
of the quadratic and the linear terms of Eq|7]with respect to {u'^)~^ . 



FIG. 6: Top: position t* of the minimum slope of MSD in log-log 
plot. Notice that for polymers t* is nearly constant , whereas for BM 
it increases with the structural relaxation time. Bottom: ratio of the 
DW factor as taken by ISF (Eq[20ll and MSD (Eq[T9ll. 



ticipates that for slowed down states, where the non-gaussian 
deviations are large (see Fig|2]i and dynamic heterogeneities 
are present, the single-monomer relaxation time Tq, and the 
collective relaxation time over a region with size ^ R^e, r^e, 
cannot be obviously related to each other. We will not analyse 
further the rotational chain dynamics and henceforth we will 
focus on the structural relaxation. 

States with different density, chain length and interac- 
tion potential are included in Fig Q corresponding to differ- 
ent degrees of anharmonicity, i.e. non-linear temperature de- 
pendence of the DW factor, and then to different fragilities 
imillllllllillillllei. The scaUng of the sti-uctural 
relaxation time in terms of Eq|7] shows that both the aver- 
age value a? and the spread cr(j2 of the square displacements 
needed to overcome the energy barriers are not affected by the 
anharmonicity. These parameters are also not affected by the 
connectivity, since the master curve collapse data of polymers 
with different chain lengths and BM. The best-fit value of the 

average is o? = 0.35, consistent with both the observation 
that (r^(t — Tq))^/^ < 0.5 (see FiglB and the well-known 
result that the atomic MSD during the structural relaxation is 
less than one atomic radius (~ 0.5 in MD units) jH]. 

The concavity of the master curve in Fig|2]is due to ~ 
0.25 7^ indicating the distribution of the displacement re- 



quired to overcome the energy barriers. Let us show that 
the concavity is a signature of the heterogeneity of the struc- 
tural relaxation. In fact, the magnitude of the ratio of the 
quadratic and the linear terms of Eq|7]with respect to {v?)^^, 
TZ = / 4:0^ (u^) , discriminates two different regimes. If 
TZ < 1 (large DW), the quadratic term is negligible and the 
displacement distribution is not observed being replaced by an 

—1/2 

effective step length , i.e. the dynamics is homogeneous. 
If 7?. > 1 (small DW), the displacement distribution shows up 
and a heterogeneous mobility distribution is anticipated. In- 
deed, on approaching the glass transition, a spatial distribution 
of mobilities develops with increasing non-gaussian features 
m [Tsl, [63ll . being characterized by the maximum a2 max of 
NGP ||63I1 . For the polymer states in Fig[Tl with NGP shown 
in Fig|2] the relation between a2 max and TZ is shown in the 
inset of FiglT] It is seen that, when TZ exceeds the unit value, 
a2 max increases exponentially. The same is observed for the 
BM states (not shown). Notably, the inset of Fig [T] reduces to 
an activated law for strong glassformers where (u^) is nearly 
proportional to T; this law has been observed for silica |63]. 



2. Other systems 

Eq|7] with the best-fit parameters from Fig|7] offers the op- 
portunity to find the DW factor (u^) at the glass transition 
of the model polymer and BM system. At the glass transi- 
tion Ta = Tag = lO^s in laboratory units |1] which cor- 
responds to Tag = 10^^ — lO^"' in dimensionless MD units 
(the time unit corresponds to 1 — IQps 1.60.1). Eq|7] yields 
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FIG. 8: Scaling of the structural relaxation time vs the reduced DW 
factor of polymers , BM (present work) and other glassformers . The 
continuous black line is Eq|23] The colored curves bound the accu- 
racy of Eq[23]and correspond to the two definitions (u^) = (r^(t = 
0.6)) (magenta, (Wg)^/^ = 0.134(1)) and (u^) = {r^{t = 1.4)) 
(orange, (y^)^^'^ = 0.122(1)X See text for details. Data taken from 
ref.(3|, IMl (Si O2), ref. Ill (icosahedral glassformer), ref. (3^ 
(Nanoparticle/polymer mixture), ref. f3'^ (OTP, data rescaled to the 
MD units of the present study). 



FIG. 9: Scaling of the structural relaxation time (in MD units) 
vs. the reduced DW factor by considering the experimental results 
concerning several glassformers. The grey area marks the glass tran- 
sition. The continuous black line is the scaling curve Eq[23] The 
numbers in parenthesis denote the fragility m. The error bounds on 
the scaling curve are the same ofFigH] Data sources are listed in 
ref.m. Ge O2 data from ref.lfsl. 



and using a, f3,j from Fig |7] yields: 



a 



(Wg)^/^ = 0.129(1). This estimate compares well with other 
related ones. First, let us consider the ratio between the vol- 
ume that is accessible to the monomer center-of-mass and the 
monomer volume is vq ^ {2{u^)^^^)^. One finds: 



va - 0.017 



(21) 



Flory and coworkers proposed that the glass transition takes 
place under iso-free volume conditions with the universal 
value vq ^ 0.025 liTsil . Furthermore, an extension of the ACS 
model (leading to the HW equation) predicts that, just as for 
a crystalline solid ll64ll , there is a Lindemann criterion for the 
stability of glasses, namely the ratio / — (u^)^/^ /d, where d 
is the average next neighbor distance of the atoms in the lat- 
tice, is a quasi-universal number (/ = 0.1) iITtIi . Our MD data 
yield 



/ 



{MD} 



0.12-0.13 



(22) 



where d is taken from the monomer radial distribution func- 
tion. /(*^^) is close to / = 0.129 for the melting of a hard 
sphere fee solid [64]. 

The knowledge of (u^) allows one to cast Eq|7]in the re- 
duced form: 



logTQ 




(23) 



-0.424(1) 
21nl0(u2) 



81nl0(u2)2 



1.62(6) 



= 12.3(1) 



(24) 
(25) 

(26) 



The above a, (3, 7 values much relies - being the largest MD 
dataset - on the evaluation of the DW factor of polymers by 
setting t* = 1 in Eq[T9] The uncertainty on t* (±0.4, see 
FigdJ leads to an error on (u^) and then on f3 and 7. One 
finds (u2)i/2 = 0.134(1) and (u2)i/2 = 0.122(1) for the two 
extremes (m^) = {r'^{t = 0.6)) a nd (m^) = {r'^(t = 1.4)) , 
respectively. By replacing in Eqs. I25l26l the extremes values 
of {Ug), the bounds setting the accuracy of Eq|23]are found. 

FigScompares Eq|23]to the results concerning the polymer 
and BM systems as well as other model glassformers. (u^) 
values of the latter were evaluated by the extrapolation tech- 
nique described above, i.e. fitting the raw data by Eq|7] and 
extrapolating the curve to Tag = 10^^ — 10^* in MD units. 
It is apparent that, within the accuracy, the scaling procedure 
works well also in these model glassformers. 



3. Experiments 

Eq|23]is well-suited for comparison with the available ex- 
perimental data. It is important to note that, even if Eql23lis 
derived in terms of {u^) msd = {u^)^ in view of the constant 
ratio {u^)isf/{u^)msd (seeFig|6ll, it also holds for {v?)isf, 
the quantity which is usually provided by the experiments. 
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FIG. 10: Correlation plot between fragility and DW factor at the 
glass transition. 



FigHjcompares the master curve Eq|23|with the experimen- 
tal data on several glassformers and polymers in a wide range 
of fragility. It covers a range of relaxation times from pi- 
coseconds to almost one year. The scaling in Fig|9| cannot 
be ascribed to (wp which weakly correlate with the fragility 
TO, see Fig[TOl Instead, it shows that both the reduced mean 
square displacement a^/ (u^) to overcome the energy barriers 
and the related spread (t„2 / (Wg) are fragility-independent, and 
then also the curvature of the master curve. 

The experimental data in Fig. [3] were collected by changing 
the temperature. In this respect, the universal scaling of Figl3] 
proves that the well-known increasing deviation of (m^(T)) 
from the linear temperature dependence of the harmonic be- 
havior by increasing the fragility index to |0, [isl 113, IH, Hit] 
just mirrors the corresponding increasing bending of Tq (T) 
vs Tg/T in the Angell plot |1] from the glass transition re- 
gion up to the liquid state. However, the glass transition may 
be reached under isothermal conditions also by increasing 
the density or the connectivity (here expressed by the chain 
length) r65'l. Our MD results highlights the correlation of 
structural relaxation and vibrational dynamics also for these 
alternative routes which awaits experimental confirmation. 



Comparison with other scaling procedures, Lindemann 
criterion 



Several scaUng and correlation plots of the structural relax- 
ation of glassforming systems were reported. They belong to 
two classes: one class, as the present approach, considers the 
relaxation times (or viscosity) of states close to and/ar from 
the glass transition |67l III llf]. We will refer to that as 
Ta-plots. The other class considers the fragility, i.e. the be- 
haviour close to the glass transition iB 1251 IItI |28i l29[ [toIi . 
These will be referred to as fragility-plots. Customarily, the 
alternative scaling procedures considered only data where the 
approach to the glass transition occurs by changing the tem- 
perature at ambient pressure, whereas the robustness of the 



Ref. 


r^"" (ns) 






Polymers 


Adjustable 
Parameters 


PW 


~ 0.001 


20 


191 


Y(6) 


1 


[66] 


~ 0.001 


32 


160 


Y(3) 


3 


[67, 68] 


~ 1 


25 


102 


N 


1 


[67, 68] 


~ 1000 


71 


174 


Y(8) 


1 


[69J 


~ 10 


53 


124 


N 


1 



TABLE 1: Comparison of the Tc-scaling procedures in literature with 
the present work (PW). rrimin and rtimax are the minimum and the 
maximum fragihty of the glass formers considered. For < r™'" 
the scaling fails. All the scaling procedures include the glass transi- 
tion region. 



Ref. 


Data 






Polymers 


Adjustable 
Parameters 


Exp 


Sim 


PW 


184 


120 


20 


191 


Y(6) 


1 


[70] 


4 




20 


90 


N 


2 


[25] 


10 




20 


87 


Y(2) 


1 


[28] 


15 




20 


100 


N 


2 


[22] 


24 




20 


160 


N 


1 



TABLE II: Comparison of the fragility-scaling procedures in litera- 
ture with the present work (PW). 



present DW scaling to pressure changes was validated by MD 
simulations. [68] 

A comparison between the present analysis and other Ta- 
plots is presented in Table U which lists the number of ad- 
justable parameters to build up the master curve, the short- 
est relaxation time t™"' below which the scaling fails, the 
fragility range being covered and the possible inclusion of 
polymers. It must be pointed out that the DW scaling adjusts 
only the conversion factor between the MD and the actual time 
units, i.e. the vertical shift factor. Within the eiTors, this factor 
is nearly independent of the system, with the notable excep- 
tion of -8203 (2d-sheet structure) |32[]. 

Fragility plots assume that fragility, i.e. the slope of the 
curve logTc vs. Tg/T at Tg, is a distinctive characteristic of 
glass-forming systems. As a consequence, they involve much 
less data than Ta-plots. A comparison between the present 
scaling and some fragility-plots is presented in Table HIl 

Recently, Niss et al considered the HW relation, EqfTj lfTTl] . 
By assuming that the intermolecular distance scales with the 
density as p^^^^, Eq[T]was recast as : 



2/3' 



(27) 



where pg and C are the density at the glass transition and a 
constant, respectively. It was concluded that, if one defines 
the glass-transition Lindemann ratio as /'^^ = p'g^^{Ug), the 
latter is system-dependent, i.e. it is not universal. It is in- 
teresting to investigate the quantity /'^-' for the polymer and 
BM models under study. It will be shown that, consistently 
with ref. lTll . /'^^•' is system-dependent. This suggests that 
/'^) is a less promising definition of the Lindemann ratio 
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FIG. 1 1 : Residues of the best-fit of log Tq vs. {u )"^) with Eq[7] 



(blue dots) and log Tq 



-2/3, 



with Eq|28](red dots). 



The ratio between the two standard deviations crjv/crpw is indicated. 
Each colored band spans ±cr. 



than f(^^) which, according to the present simulation, de- 
pends very weakly on the system, see Eq|22] It must be re- 
minded that our MD data set correspond to different kind of 
polymeric and BM systems in that different interacting poten- 
tials, Eq[TT] are considered and, for polymers, different chain 
lengths. To test the quantity /(^'' , the correlation plot between 
logTa and was fitted with: 



logTa 



{N) 



Ci + C2p-^'HuY' + C^p-^r^u')-'' (28) 



/«2x 



The above form is analogous to Eq|2]and suitably general- 
izes Eq|22]to account for the bending of the plot of log Tq vs. 



-2/3/ 



^ due, in turn, to the bending of log vs 



(see Fig|7j. Eq|28] assumes that the density scaling of the 

-^1/2 

characteristic length scales fulfills the ansatz , (T(j2 



1/2 



(X 

p"^!^ . Instead, Eq|7]takes both quantities as constant. FigfTT] 
compares the residues of the fit of log vs with Eq|7] 

(see Fig|2ll and the fit of logT„ vs ^^2/3^^2^-1 ^-^^ EqESl 
Both fits have the same number of adjustable parameters. It 
is seen that the residues are structureless, i.e. the fits have 
equal accuracy. However, the discrepancies from Eq|28] ex- 
hibit standard deviation which is larger than the one, apw, 
of the deviations from Eq|2] Note also that the deviations 
from Eq|28] increase with r^. When an extrapolation proce- 
dure analogous to the one outlined in Sec llVB 2l is followed 



when logTa is plotted vs. the quantity p '^^'^{v?) ^ results in 
a /(^) value with larger uncertainty, i.e. less "universal", than 

j(M_D) 



V. CONCLUSIONS 

The paper presents a thorough analysis of the scaling be- 
tween the long-time relaxation and the caged dynamics. MD 
simulations of both a model polymer system and a binary mix- 
tures were carried out by varying the temperature, the density, 
the potential and the polymer length to consider the structural 
relaxation as well as both the rotational and the translation dif- 
fusion. They showed the existence of different physical states 
exhibiting coinciding transport and relaxation from the end of 
the ballistic regime through the diffusive one. This points to 
a link between the short- and the long-time dynamics which 
is evidenced by the master curves found by correlating the 
DW factor with the structural relaxation time and the chain 
rotational diffusion. An analytic model of master curve 
is developed, leading to Eq|7] which fits nicely with the MD 
results on polymers and BM. Notably, the model does not pre- 
dict the existence of physical states yielding the divergence of 
Tot- By using suitable reduced units, the MD scaling on 
polymers and BM was extended to include MD data on other 
systems as well as the experimental data on several glassform- 
ers and polymers in a wide range of fragility by covering a 
range of relaxation times from picoseconds to several days. 
The scaling in terms of the DW factor compare favourably 
with other scaling procedures. In particular, the density scal- 
ing of the characteristic length scales according to the ansatz 

cr„2 cx p^^l^ is not supported by the present simula- 



>l/2 



to derive p 



,2/3/ 



,.e. /(^), the poorer collapse of the data 



tions. The study suggests that the equilibrium and the moder- 
ately supercooled states of the glassformers possess key infor- 
mation on the huge slowing-down of their relaxation close to 
the glass transition which, according to our simulations, ex- 
hibits features shared with the Lindemann melting criterion 
and the free-volume model. 
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